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Abstract 



Simulated minimum cross-section histograms of breaking Aluminum 
nanocontacts are produced using molecular dynamics. The results allow a 
new interpretation of the controverted conductance histogram peaks based on 
preferential geometrical arrangments of nanocontact necks. As temperature 
increases, lower conductance peaks decrease in favor of broader and higher 
conductance structures. This reveals the existence of supershell structures 
favored by the increased mobility of Al atoms. 

PACS numbers: 73.40. Jn,61. 16. Ch,73.23.Ad 
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Electron transport in metallic nanowires having cross-sections formed by few atoms ex- 
hibit a rich quantum phenomenology since their characteristic sizes are of the order of the 
electron Fermi wavelength (Ap). Transverse confinement of electrons in the nanocontact 
region only allows few propagating modes N, giving rise to quantization of the conductance 
(assuming perfect transmittance per channel) G = N x Go, where Go = 2e 2 /h is the con- 
ductance quantum Although conductance quantization (CQ) was first seen in 2DEG 
devices ||, the study of this phenomenon in metallic nanocontacts has received much at- 
tention because of its technological implications in mesoscopic swithching elements and new 
nanoelectronics devices |JJ. 

Metallic nanocontacts have been built with scanning tunneling microscopy (STM) [f|||, 
with the Mechanically Controllable Break Junction (MCBJ) method |6|]7}] as well with two 
plain macroscopic wires M. Nanocontact formation and later breaking are usually carried 
out in a controlled way using piezoelectric actuators. The conductance evolution is followed 
by applying a bias voltage to the nanocontact. In general, conductance jumps appear during 
the nanowire breaking process and correspond to changes of the transmission probability of 
the propagating modes or sudden variations of N. 

Since every single contact breaking experiment has its own conductance evolution, the 
study of CQ has been usually addressed from a statistical point of view, elaborating con- 
ductance histograms (CH) to get information on the electrical behavior of nanocontacts. 
In general, CH show a well-defined peaked structure, with higher conductance probabilities 
close to integer values of Go- This trend has been usually interpreted as the evidence of the 
CQ in different metallic species . CH have been used to evaluate changes in transport 

properties when modifying external parameters such as the magnetic field ||, the chemical 
environment [10], temperature [11,I2]1, or the applied bias voltage [13 . 



Although CH has become a standard tool to study the electronic transport in metallic 
nanocontacts, there is no solid theoretical background justifying their use. Any realistic 
interpretation of CH should take into account that there is a strong coupling between elec- 
tronic and mechanical properties 
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since conductance and force jumps are correlated, 
through crystalline rearrangements inside the nanocontact as predicted by Molecular Dy- 
namics (MD) simulations ]TS|. Therefore, CH give information of the conductance, while 



the measured histogram weights are related to the stability of a given conductance situation. 
On the other hand, while for monovalent species an interpretation of CH in terms of number 
of atoms at the neck seems straightforward [16[], for species with higher chemical valence 



the scenario is less clear. Scheer et al. WR showed that three modes contribute to electron 
transport in aluminum monoatomic contacts whose conductance is close to Go- Similar 
results for Pb and Nb demonstrate the existence of a correspondence between valence and 
the number of conduction channels participating in one-atom contacts |T7| , |T8|1 . For Al, CH 
constructed from many MCBJ experiments fl9j clearly show peaks close to G/Go ~0.8, ~ 
1.8, ~3.0, and ~4.0 at T=4K. A quite similar experimental series of peaks for Al has been 
reported recently [p0| . These results have been used to reconsider the role of CH to indicate 
CQ features. In particular, it has been suggested that peaks in aluminum CH originate in 



favorable atomic configurations appearing during nanocontact stretching fT9| , |20 



In the past, within the framework of the free-electron model together with conductance 
calculations, theoretical CH have been obtained by considering a simplified ad-hoc nanowire 
dynamics [J21]J. A realistic description of the dynamics of the metallic nanowires breakage 
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can be obtained by means of first-principles calculations |2~2fl. Other authors used classi 



cal Molecular Dynamics accompanied with conductance calculations [p3~25l- However, to 
construct CH, both approaches are impracticable because of the prohibitive computational 
cost. 

In this letter, we have adopted a different perspective. Our goal is to construct atomic 
configuration histograms of the neck instead of conductance histograms. In this way we 
will be able to know whether there are favorable atomic configurations behind the peaked 
structure obtained in CH. We have studied the atomic configuration histograms of aluminum, 
a trivalent metal, which shows a well defined and reproducible peaked structure [iT5,2TJ. 



For the statistical study of breaking nanocontacts we performed MD simulations using the 
embedded atom method (EAM) |26| as in previous works [|i~5p 4"f. In particular, we have 
used state-of-the-art many-body EAM interatomic potentials for aluminum able to fit bulk 



and surface properties p7| . Simulations were carried out in a wide range of temperatures 
(4K-450K) in order to assess the influence of temperature T on the configuration histograms. 
The temperature is controlled during the simulation using conventional scaling of velocities. 

As starting point of our MD simulation we consider a parallelepiped supercell containing 
1008 Al atoms ordered following a fee crystallographic structure. The lattice constant is 
taken to be equal to the experimentally measured value (4.05 A). Aluminum atoms were 
initially distributed in 18 layers perpendicular to the (111) direction containing 56 atoms 
each [see Fig. la (left) for illustration]. The direction (111) (z axis in our simulation) corre- 
sponds to that in which the contact will be elongated until breaking, although we also have 
considered other orientations discussed below. Each single MD simulation is carried out fol- 
lowing three stages. First, we relax the initial bulk-like ordered system for 2000 iterations, 
imposing periodic boundary conditions (PBC) in the x and y directions of the sample. The 
time interval per iteration step in all our simulations is dt = 10~ 14 sec. At the end of this 
equilibration step the system undergoes a relaxation characterized by a contraction along 
the z axis. In the final thermally equilibrated structure two bi-layer slabs are defined at the 
top and bottom of the supercell. All atomic positions inside these slabs are frozen during 
subsequent MD stages, defining the bulk support of the nanocontact during the breaking 
process. In a second stage, the system is again thermally equilibrated during 3000 itera- 
tions, but now PBC are preserved only in the x and y directions of the frozen slabs. The 
other atoms, defining a free nanowire evolve to reach a new equilibrium configuration. For 
temperatures higher than 450K, we found that we are not able to stabilize the nanowire 
geometry, which spontaneously breaks. This is consistent with the decrease of the melting 



temperature noted for metallic nanowires as their radius decreases [28]. The third stage cor 



responds to the nanowire stretching, separating both supporting frozen slabs at a constant 
velocity of 2 x 10 -4 A per iteration step. The full determination of atomic positions during 
the elongation process allows the determination of an effective surface for the nanocontact 
and a minimum cross-section Sn time evolution The quantity Sn is computed as in ref. 
f2~5fl (in units of numbers of atoms) and is approximately equal to the Sharvin conductance 
of the nanocontact in units of Gq. In order to reduce the computational effort, we evaluate 
Sn each 10 MD iterations. 

An illustration of the nanocontact evolution during elongation Ad is depicted in Fig. 
la. Fig. lb shows the evolution of Sn as a function of Ad for four different realizations 
(i.e. each relaxed pallelepiped gives rise to a different nanowire evolution until breaking 
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takes place). The jumps in Sn curves are correlated with the jumps of the resulting force 
on the frozen slabs as observed previously @,[25[]. By accumulating conductance traces we 
can construct minimum cross-section histograms H(Sn)- The number of single nanowire 
breaking realizations N s considered is such that the obtained histogram peak structures are 
unchanged under further averaging (see Fig. lc). In our calculations we have found that 
N s =100 provides a reliable H(Sn) curve. 

In Fig. 2 we plot the histogram H(Sn) for three different temperatures. For low tem- 
peratures (T = 4K) it is clear that H(Sn) presents a well-defined peaked structure with 
maxima located at Sn = 1,2,3,4, etc. Therefore, the statistics of the nanocontact at the 
narrowest section exposes favorable atomic configurations during its stretching. Although 
our histograms do not contain direct information on conductance, the resemblance with 
experimental histograms fl9lJ2"0|| is remarkable. Assuming that one Al atom provides a con- 
ductance close to Go (including components from three different channels) the interpretation 
of H(Sn) as a conductance histogram is close to the observed results. 

For increasing temperatures, we found that the histogram dramatically changes loosing 
its peaked features for low Sn values, whereas a different structure appears for high Sn 
values. This trend is similar to that found in Na conductance histograms reflecting the 
appearance of supershell high stability structures |T2[. In our case, the temperature needed 



to find such supershell configurations is higher than in Na wires due the higher cohesion 
energy of Al (1.13 and 3.36 eV in Na and Al bulk, respectively). Although the explored 
region only corresponds to Sn < 30 the presence of peaks at Sn — 5, 8, 12, 17, 24 is evident. 
The shell structure becomes observable at higher temperatures because of the increased 
mobility of the atoms. This allows the system to explore many atomic configurations in 
order to efficiently find a suitable local energy minimum. 

By inspecting Fig. 2 (where H(Sn) is shown in a range similar to that used in experimen- 
tal situations), it is clear that temperature has a dramatic effect on the favorable geometrical 
configurations. The higher the temperature the larger the background structure, giving rise 
to a histogram with a less defined peak structure. In particular it is very striking that the 
peak at Sn = 1 is still well defined at T = 300K while the peak Sn = 2 decreases its relative 
weight with respect those found in Sn = 3 and 4, in comparison to the 4K case. At T=450K 
the Sn = 1 peak almost vanishes and the broad structure around Sn = 5 appears, whereas 
intermediate peaks at Sn = 2,3 and 4 disappear. Although these predictions correspond to 
cross-section histograms, we believe a similar behavior is to be expected in true CH. 

In actual experiments the nanowire orientation is not well defined, and the histogram 
would contain components coming from different nanowire orientations. In order to analyze 
these effects we have studied the histograms H(Sn) at T = 4K for two other orientations 
of the stretching direction: (Oil) and (001). These results are depicted in Fig. 3. It is 
worth noticing that the histograms have now a more complicated structure but also reflect 
the presence of favorable atomic configurations. We can then expect that in experiments 
the histograms can be understood in terms of a combination of the histograms appearing 
in Fig. 3. However it is clear that a structure with peaks at Sn = 1 and 2 should be 
present. For Sn > 3 the peaked structure should become less clear and would depend 
on the experimental conditions (as it happens when comparing two different experimental 
conductance histograms in Al [pl||||]). Another relevant observation is that for less stable 
configurations, as those considered when stretching along (011) and (001) directions, the 
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shell structure tends to appear at lower temperatures [see, for instance, the peak at Sn = 5 
formed in the (Oil) case]. Thus, orientation effects will also contribute in bringing out 
supershell effects. 

In conclusion, we have obtained, for the first time, minimum cross-section histograms for 
Al nanowires breaking MD simulations using state-of-the-art EAM potentials to describe 
Al-Al interactions. At low temperatures we found that there exist favorable atomic ar- 
rangements as suggested some time ago by experimentalist []I9|,Z0|. This provides a new 



way to interpret experimental conductance histograms on the basis of MD simulations. For 
increasing temperatures, the system is able to explore more atomic configurations where low 
conductance regions show less structure and appearing well defined high stability regions for 



higher conductance values appear, as reported for Na |12| . In the low conductance region 



it is found that at T = 300K the histogram peak at Sn = 2 decreases its relative weight 
with respect the peak Sn = 3. That is, at 4K the peak height sequence is 1, 2, 3, 4 whereas 
at 300K the sequence change to 1, 4, 3, 2. For T=450K the peak at Sn = 1 almost van- 
ishes, and the presence of supershell structure becomes evident. Finally, for other nanowire 
stretching directions in the low temperature regime we also found that there exist favor- 
able configurations, although the histogram presents a more complex structure. However 
peaks at Sn = 1 and 2 remain. Nanowires with lower stability tend to form the supershell 
structure at lower temperatures than those with high initial stability. 

We acknowledge J.J. Saenz, P. Garcia-Mochales, C. Urbina, R. Paredes, J.R. Villarroel, 
and A. Garci'a-Martm, for helpful discussions. This work has been partially supported by 
the CSIC-IVIC researchers exchange program and the Spanish DGICyT (MEC) through 
Project PB98-0464. 
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FIGURES 



FIG. 1. (a) Typical nanocontact configurations during the elongation of the sample, (b) Mini- 
mum cross-section Sn evolution for four different Al nanowires at T = 4K. The stretching direction 
is parallel to the (111) direction, and the arrows indicate the corresponding Sn for each configura- 
tion, (c) Minimum cross-section histograms H(Sn)- Dark gray, gray and black curves correspond 
to histograms obtained with different number of realizations (35, 70 and 100, respectively). 

FIG. 2. Minimum cross-section histograms H(Sn) for stretching of Al nanowires along the 
(111) direction at T=4K, 300K, and 450K. The right panel is a zoom of the range < Sn < 8. 
The curves are normalized to the number of single breaking events (100 per temperature value). 

FIG. 3. Comparison of the minimum cross-section histograms H(Sn) of aluminum nanowires 
at T = 4K stretching along three different orientations (111), (011) and (001). The last two 
histograms are constructed with 70 breaking events and nanocontacts contain 1000 and 980 atoms, 
respectively. All histograms are normalized to the number of realizations 
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